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Abstract The quantum entanglement is an important feature of many pro¬ 
tocols in the field of quantum computing. In this paper we evaluate a level of 
entanglement in short qutrit chains. This evaluation is carried out with use 
of the CCNR criterion and the concurrence measure. We also present some 
explicit formulae describing the values of CCNR criterion and concurrence for 
exemplary short spin chains. Utilizing the obtained results, we indicate that 
analyzing the level of entanglement allows to detect the noise or deviation in 
the transfer process, in comparison to the perfect transfer where only operation 
realizing transfer is present. 

Keywords quantum information transfer • qubit/qutrit chains • entangle¬ 
ment • ccnr criterion • concurrence • correctness of transfer • perfect transfer • 
numerical simulations 


1 Introduction 

A phenomenon of quantum entanglement |12j is a fundamental feature of quan¬ 
tum systems. A presence of quantum entanglement in an examined system 
indicates system’s quantum nature. The entanglement is also present in spin 
chains which became one of the essential elements of quantum physics and 
quantum computing m. m, ca¬ 
using spin chains to transfer an information was shown for the first time 
in [^. An issue of the perfect transfer was widely discussed in 0, m, m- 
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An important problem of qubit and qudit states’ transmission was raised in 
n, 0, m- The next issues are: creating, detecting and analyzing the level of 
entanglement in two adjacent chain’s nodes m. la. The questions connected 
with the presence of entanglement in spin chains are still an active research 
field. In [57], |5S], [55], [30] the analysis of entanglement in spin chains, also 
for qutrit chains, was presented. 

The main function of spin chains is a transfer of an information. A level 
of entanglement in a spin chain changes during the transmission of quantum 
state/information. In spite of difficulties associated with formulating the en¬ 
tanglement levels, it is possible to evaluate the entanglement e.g. in spin chains. 
In this paper the level of entanglement is estimated for short qutrit spin chains. 
The presented results may be used to verify if the process of transfer is correct, 
because it seems that the level of entanglement can be treated as an invariant 
for the transfer protocol in a spin chain. 

The paper contains the following information: in Sect. (0) a form of Hamil¬ 
tonian for performing a XY-like transfer protocol for qutrit and qudit chains 
is presented (including a short note concerning properties of the perfect trans¬ 
fer). There is also an algorithm describing the realization of transfer protocol 
with a 7 condition playing a role of the invariant expressing a level of entan¬ 
glement during the transfer. In Sect. (|3|) a chosen criterion for entanglement 
detection and the concurrence measure are presented. 

Sect. (0) contains the results of the experiments for detecting entanglement 
and calculating the values of concurrence measure. 

In Sect. m we show the influence of noise on the data transfer process. 
There is also an average value of Fidelity computed for a channel where dis¬ 
tortions are presented by a phase-damping method. 

A summary and conclusions are presented in Sect. (ED- 


2 Hamiltonian for a Transfer Protocol in a Qudit Spin Chain 

In this section we define a Hamiltonian which will be used to realize the 

perfect transfer of quantum information in qutrit and qudit chains [Sj, m, 
0 for entanglement creation between chosen points of a chain. The form of 
Hamiltonian, given below, is naturally suitable for transfers discussed in this 
work - i.e. for transmission of information in qutrit spin chains. 

The SU{d) generator is utilized to create the XY-like Hamiltonian for 
qudits. The construction of suitable SU(d) generator is given in Appendix A. 
Assuming that each qudit has the same freedom level d > 2 (the qudit is 
defined in a similar way to qubit, however a computational base for qudits is 
expressed with d orthonormal vectors - in case of qubits, the base contains 
two orthonormal vectors): 



( 1 ) 


(2,2+1) G-C(G) 
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where Ji is defined as follows: J, = for 1 < fc < j < d and 

are SU{d) group operators defined by (15^ applied to the (*)-th and (i + l)-th 

qudit. The Hamiltonian © will be also called the transfer Hamiltonian. 

The state transfers, studied in 0 , 0 , m use Hamiltonian H which have 
the following property 

N 

[iJ,^%]=0, (2) 

i=l 

where Z represents the sign gate for qubits. This means that spins are pre¬ 
served and dynamics generated by H is divided into series of subspaces denoted 
by the number of qubit in state |1) - see in [T3]. In the case discussed here, it 
is not hard to show that 

N 

= ( 3 ) 

i=l 


for 1 < r < (d — 1), so the equation (|3l) generalizes the situation mentioned in 
the equation ([ 2 ]) - preserving spins and separating dynamics into subspaces. 

It is necessary to add that an appropriate unitary operator for transfer 
operation is determined by the equation 


Uf = e 




(4) 


where t represents evolution time and i represents imaginary unity. The symbol 
U represents the unitary operation which performs the transfer protocol. 

A very important issue raising in the context of information’s transfer in 
spin chains is a problem of the perfect transfer. If the symbol 1 denotes the 
initial node of the chain and N stands for the final node then a transfer is the 
perfect transfer when for the time t we have: 

|(iV|e-^""^‘|l)| = I. (5) 


In such a situation there were no distortions and the final and initial states 
are the same in the meaning of Fidelity measure. 

Introducing a definition of the transfer protocol based on a Hamiltonian 
and a unitary operator allows to describe the transfer (Fig. [T]) as an algorithm 
(or structural quantum program [23], [I6j'l. A very important issue is the use 
of 7 condition as an invariant for the protocol. The invariant 7 is based on a 
function calculating the level of entanglement. 

Let E{\<I'{t))) be the function to determine a level of entanglement for a 
chain under perfect state transfer. The CCNR criterion, negativity and con¬ 
currence seem to be good candidates to play this role. 


Remark 1 Calculating the level of entanglement in multiple qubit, qutrit and, 
particularly, qudit systems is a problem which is still unsolved for mixed quan¬ 
tum states (however it is possible to give lower bound for concurrence for 
mixed states 0)- Due to this fact, in the notation of invariant 7 for perfect 
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t ; time variable 

> 0 ; number of step 

j Hamiltonian for path L with I vertices) 
f/t = ; unitary operator 

\'ip{to)) ; initial state of chain 

i : = 0 ; 

while i < N do 
begin 

|f/j(ti+l)> := Utl^jiU)) 
i : = i + 1 

{y ■. E{mi))) ^ Eiint)})} Ai < N 

end ; 

h ■■ E{mi))) E{\Em Ai > N} 

Fig. 1 The data transfer process in a spin chain showed as an algorithm expressing the use 
of a unitary operator on a quantum register. The condition 7 plays the role of an invariant 
and describes the level of entanglement in the register |' 0 (ti)) where ti represents the time 
in a step number i 


l^o) IV'l) IV’2> |V'ti-2> llpn-l) I'tpn) 

Fig. 2 The picture of quantum information transfer’s realization in a spin chain for a single 
unknown qudit state |'0) in the path of length 1. The state represents the initial state 
and state lipn) denotes the final state obtained after execution of transfer protocol 


state protocol (algorithm) - 7 : E{\'tlj{ti))) = E{\<E{t))) - sign = expresses 
that the level of entanglement is comparable with used measure (e.g. CCNR 
criterion or concurrence). 

In this paper the transfer protocol is used to realize the transmission of 
information through the path which length equals N. In general the process 
of transfer involves an unknown qudit state (where d is a freedom level): 

IV') = Q;o| 0) + ai|l) + ... + ad-ijd — 1) where G C and *< d. (6) 

However, in our examples described in further part of this work, we transfer 
only the qutrit state, i.e. d = 3. 

Fig. [2] shows the scheme of transfer protocol’s realization. The whole pro¬ 
cess is divided into discrete steps which are realized by operator Ut- 


3 The CCNR Criterion and Concurrence 

In this part a basic information about CCNR criterion and concurrence will be 
recalled. The mentioned methods may be used as measures to estimate a level 
of entanglement in a bipartite system. The CCNR criterion uses two auxiliary 
operations called vectorization and realignment. 
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(a) 


(b) 


IV’in) 





Fig. 3 The examples of quantum circuits realizing the transfer protocol for N nodes. In 
the case (a) we assume that the whole process will be realized by one unitary operator T, so 
this operation needs an access to all chain’s nodes. The case (b) shows the quantum circuit 
where the transfer process is realized by local operators Ti correlated with two adjoining 
chain’s nodes 


The vectorization of a given matrix A = [aij]mxn is represented by the 
operation vec : Amxn —t C™", and 

'Uec(^) — [^11? ■ ■ ■ 1 ^ml-j ^12; • ■ ■ ; ^m2; ■ ■ ■ 5 ^ln^ • ■ • ; ^mn] (7) 

The realignment is based on the vectorization of matrix and applies to the 
scalar product of two matrices: 

R{A®B)=vec{A)^ ■vec{B). (8) 


The realignment of matrix’s elements - the initial matrix dimensions are 2x2 
- in a space H = Ha ® Hb may be illustrated by the example: 


/ 011 012 

Ol 3 014 \ 


/ Oil O21 

O12 022\ 

021 022 

023 024 

II 

O31 O41 

O32 O42 

031 032 

O33 O34 

Ol3 O23 

Ol 4 024 

\ 04 i 042 

O43 O44 / 


\033 O43 

O34 O44 / 


The operation of realignment is realized similarly for larger matrices, e.g. 
for matrices 3x3: 



/ Oil 

O21 

O31 

O12 

O22 

O32 

O13 

O23 

O33 


041 

O51 

061 

O42 

O52 

062 

O43 

O53 

063 


071 

081 

O91 

O72 

082 

O92 

O73 

083 

O93 


014 

024 

O34 

Ol 5 

O25 

O35 

O16 

026 

036 

R{p) = 

O44 

O54 

064 

O45 

O55 

065 

046 

056 

066 


O74 

084 

O94 

O75 

085 

O95 

076 

086 

096 


Ol 7 

O27 

O37 

O18 

028 

038 

O19 

O29 

O39 


O47 

O57 

067 

048 

O58 

068 

O49 

O59 

069 


\077 

087 

O97 

078 

088 

098 

O79 

089 

O99 


The realignment allows to specify the criterion which is able to detect and 
track the entanglement’s behaviour during the transfer process. This method 
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is called CCNR criterion m- Generally, the CCNR criterion is defined by the 
fact that if the matrix pab of a bipartite m x n system is separable, then: 

II4bII<i, (11) 

in the case Hp^sll > 1 state pab is entangled. Of course p^ stands for the 
completion of realignment operation on the state p. 

Remark 2 Generally, the CCNR criterion is not sufficient to detect all the cases 
of entanglement in any quantum system. However, it is sufficient to correctly 
detect the entanglement in qudit, and of course qutrit, spin chains. 

It should be also pointed out that the value of CCNR criterion may be 
calculated using Singular Value Decomposition (SVD): 

IIP^bII = ^ijPAB) (12) 

i=l 

where at represents a singular value of Pab s-nd Q = min(m^,n^). 

The correctness of CCNR criterion may be amplified by using the states 
of subsystems what was shown in EH: 

WiPAB - Pa<S)Pb)^\\ < Y^(l - Tr(p^))(l - Tr(p|)), (13) 

If the above condition is fulfilled, it means that the state pab is separable 
(otherwise the state is entangled). 

The concurrence measure is used to calculate a level of entanglement, as 
well. Concurrence [5] for pure state {iPab) in bipartite system d® d is defined 
as 

Cili^AB)) = ^2(1-Tr(pi)). (14) 


4 Tracing of Entanglement 

In this section we present the results of calculating a level of entanglement 
with use of CCNR criterion and concurrence. Fig.|3]shows the values of CCNR 
criteria calculated according to the equations (HU and m- The process of 
transfer was realized for the state |+), which for qubits is |+)2 state and for 
qutrits is expressed as |+)3 state: 

|+)2 = ^(|0) + |1)) or |+)3 = ^(|0) + |1) + |2)). (15) 

The value of CCNR criterion was calculated according to the equation (fT^ . 
In case of amplified CCNR criterion, the equation (TT^ was transformed to: 


\\{PAB -PA® Pb)^\\ - (1 - Tt {p\)){l - Tt (pD) > 0 


(16) 
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Fig. 4 The entanglement detection with use of CCNR criterion and amplified CCNR cri¬ 
terion for short qubit and qutrit chains. The presented values were calculated numerically 
and the transferred state was |-|-} state 


If the above inequality is true, then the analysed state is entangled. 

Naturally, during the transfer process, the level of entanglement initially 
increases and then decreases. All analysed short spin chains were divided into 
two subsystems: A and B. 

Remark 3 It should be pointed out that the subsystems A and B are treated 
as systems with a greater number of dimensions. However, it is still possible 
to detect the entanglement with use of CCNR criterion and to evaluate the 
level of entanglement with concurrence. 

Fig.[S]shows the values of concurrence measure for some exemplary chains. 
It can be observed that the entanglement’s dynamics is similar to the results 
from Fig. m 

Using the equation da allows to present the explicit formulae for calcu¬ 
lating the level of entanglement in spin chains. If an unknown qubit state: 


\ijj) = alO) -h/3|l), 


(17) 
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Fig. 5 The values of concurrence measure for short qubit and qutrit chains. The presented 
results were calculated numerically and the transferred state was |+) state 


is transferred through the short path {1 = 2), then the level of entanglement 
in this 2-qubit chain can be expressed as: 

i ( 4 q;^ + 3 / 3 ^ + 8 q:^/ 3 ^ cos( 2 a) -I- ( 3 '^ cos( 4 a)) ( 18 ) 

where a represents the duration of the transfer process. 

For a 2-qutrit chain the transfer of an unknown state: 

|V')=a|0)+/3|l)+7|2), (19) 

can be precisely described by the concurrence: 

Cd=|(a) = i (4a‘^ + 3(/3^ -f 7 ^)^ -f 8 a^ (/3^ -f 7 ^) cos( 2 a) -f -f 72 )%os( 4 a)) 

( 20 ) 

The dynamics of changes for concurrence values depends on the values 
cos(2a) and cos(4a) - both in qubit and qutrit chains. The constants preced¬ 
ing cos(2a) and cos(4a) correspond to the probability amplitudes. A sum of 
probability amplitudes’ values equals one what is easy to proof when for some 
a values: cos(2a) = cos(4a) = 1. Furthermore the values of amplitudes: 


{a^+P^)^ = l, {a^ + p^+^Y = l 


(21) 
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Fig. 6 The values of concurrence measure, for short spin chains, obtained with use of the 
formulae di and d 


both for qubit and qutrit chains of any length, because of a trace calculation 
in m- An exemplary value of concurrence measure for 4-qutrit chain is: 

Cliztio,) = ^ (^0 + (ci + C2) cos 2a + (03 + C4 + C5) cos4a+ (22) 
(c6 + C7) cos 6a + (cs + cg + cio) cos 8a + (cn + C12 + C13) cos 12a) 

where Ci = 1 and the Ci coefficients are composed of the probability am¬ 
plitudes’ values of transferred state. 

Fig. ini presents the values of concurrence measure for explicit formulae. 
The changes in the level of entanglement, illustrated as local minima, were not 
detectable with numerical calculations (see Fig. |T|). These minima correspond 
to the different levels of entanglement between separate chain’s parts, e.g. 
between adjacent qubits/qutrits. 

Remark It should be noted that usually there is no need to analyze the 
global level of entanglement for the whole chain. The character of a transfer 
in a one-dimensional spin chain causes that the entire transferred information 
will be placed in the last node in the last step of the process. It means that 
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the levels of entanglement can be calculated only in the first and the last node 
to check if the transfer process is correct. The needed quantum state, for the 
spin chain with the number of nodes greater than two, we obtain by use of 
the partial trace operation. If p represents the state of the whole chain with 
N nodes then the state in the first (denoted as 1) and the last (described as 
N) node we calculate as: 

Pin = Ti'(2,3,...,i-fc,...,Af-2,Ar-i)(p) (23) 

where Tr(.)(p) is the partial trace operation. 


5 Transfer with Noise Presence 

Naturally, the presence of noise in a transfer process causes a distortion of 
quantum state what can be verified, i.a., with use of Fidelity measure. The 
presence of noise may be analyzed in the context of Grover’s algorithm m and 
usually the distortions severely hamper the realization of a quantum algorithm 
or a quantum protocol. 

However, calculating the value of Fidelity needs the information about the 
quantum state before and after the transfer. We know that the initial state is 
separable and the phenomenon of entanglement should only be present during 
the transfer process (if the transfer is perfect, the final state is also separable). 
This issue should make us reflect that calculating the level of entanglement 
for the final state, with use one of known measures, may tell us if the transfer 
was perfect. 

The concurrence is relatively easy to compute so we will use this measure 
(HU to detect the presence of noise. 

The distortions will be presented as a phase-damping in a qutrit spin 
chains, expressed as a quantum channel in a Kraus representation. Phase¬ 
damping will be applied to the whole qutrit chain. However, for qutrits, and 
for qudits in general, the phase-damping operation does not have a unique 
representation. The model discussed in the works m, [I] is one of the 
examples of this operation. In our work we use the following operators for 
phase-damping: 


d-l 

£{p) = '£e,pe}, 


2=0 


E^ = 


i(E)Mwr 


(24) 


where the damping strength is as follows 0 < p < I and (f ^ represents the 
binomial theorem. 


Remark 5 It should be noted that expression (l24l) can be regarded as a special 
case of Weyl’s channel [5] : 

d-l 

£{p)= ^m.n(^”X'")p(X'"Z'‘)^ 

m,n—0 


(25) 
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(a) 


(b) 






Ti 


£i 






Tn-1 


£n-i 


{‘4’out) 


(c) 




Fig. 7 The diagrams of quantum circuits realizing the transfer of quantum state. In a case 
(a) the transfer is realized with use of one T gate (this operation is performed in every 
node) and the distortion is presented as an S gate. A case (b) expresses the situation where 
a transfer and a distortion are realized locally (effecting only adjacent nodes). In a case (c) 
the operators T and S are decomposed to a sequence of Ti and Si operations which are 
performed alternately 


where elements of the matrix tt satisfy the following conditions: 0 < < 1 

and Yltn n=o = 1- The operators Z and X are generalised Pauli matrices 
for the sign changing and negation operations on qudits. 


The definitions of qudit gates X and Z are: 


^|j) = |j®l), Z\j)=uj^\j), 
where w stands for fc-th root of unity: 


/2fc7r\ /2fe7r\ 2 ,rik 

= cos —— + i sm —— = e 


\ d J 


\ d 


k = 0,1,2,1, 


(26) 


(27) 


d expresses the degree of root and i denotes the imaginary unit. It is also 
assumed that 

d 

j © a = (j + a) mod d, j G N. (28) 

The distortion operator will be utilized after the operation of data transfer. 
The transfer operator is a global operator which affects all the elements in a 
short spin chain. Three basic ways to introduce the noise to a spin chain are 
presented at Fig. ©■ 

Fig. [5] depicts the values of concurrence measure for a few experiments with 
short spin chains. The diagram shows that the presence of noise significantly 
changes the level of entanglement during the transfer process (the level of 
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Fig. 8 The value of concurrence measure during the data transfer in two and more qutrit 
chains. If the noise is present, it is generated as a phase-damping with parameter p = 0.85 


entanglement should decrease in the final part of the process). The distortions 
prevent the drop of entanglement’s level, so the presence of entanglement in 
the final state means that the transfer was not completed successfully. 

We can confirm this result also by calculating the average value of Fidelity 
for the perfect transfer and transfer operators in e.g. 2-qutrit chains. The works 
m, m, m show that the average value of Fidelity may be computed as: 

^avg 

where Mk = Ut^Ek and £{p) = ^kpEj, represent the quantum channel - in 
this case these are phase-damping operators. The unitary operator Ut realizes 
the perfect transfer. The average value of Fidelity, calculated analytically, for 
a 2-qutrit spin chain; 


Fayg = ^{Sp^+y-l\+ip+S). (30) 

The average value of Fidelity equals 2/3 what is confirmed by our numerical 
experiment (see Fig. ([5])) - the average value of Fidelity, calculated numerically, 
is Fayg « 0.62702. 

At Fig. (O the level of entanglement is also presented. Its value confirms 
that the transfer process did not end correctly. 
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Fig. 9 The values of Fidelity and Concurrence measures for the 2-qutrit chain where the 
transfer process was divided into sixteen steps and the distortions were caused by phase¬ 
damping (with different parameters) 


6 Conclusions 

Calculating the level of entanglement with use of CCNR criterion and, es¬ 
pecially, with concurrence proved to be a good solution for tracking the in¬ 
formation transfer in qutrit spin chains. The concurrence measure indicates 
unambiguously the details of transfer process and due to this fact, it can be 
used as an invariant in the algorithmic description of transfer protocol for an 
unknown qubit state and also for an unknown qutrit state (generalization for 
qudit states is also available). 

Naturally, the presented analysis concerns only short spin chains and it 
should be extended for chains of any length. However, this extension is in fact 
only estimating the upper bound for concurrence value in a spin chain. The 
numerical experiments seem to imply that the length I = 2 points out the 
upper bound for concurrence in a spin chain given by the Hamiltonian CD- 

Introducing the noise in a form of phase-dumping causes a higher level of 
entanglement and a lower value of the Fidelity measure. Tracking the entan¬ 
glement level in chains with and without noise allows to study the distortions. 
We suppose that analyzing the level of entanglement in chains with and with¬ 
out noise, for example to compare the values of Fidelity measure, would be an 
interesting research issue. 

It is important to add that the concurrence measure - apart from its ana¬ 
lytical functions, e.g. to track the level of entanglement, as it was shown in the 
paper - is at present also used in the experiments [5] where the entanglement 
is produced in a two-photon system. In |25j the concurrence was utilized for 
spin-1/2 chain which can be realized physically with use of Benzene molecules 
. Therefore, the functions described in the chapter may be applied as a tool 
to verify the concurrence in the physical experiments where e.g. qutrits were 
used. 
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Appendix A Construction of Lie algebra’s generator 


In a following definition of the XY-like Hamiltonian for qutrits’ and, generally, 
qudits’ chain, the Lie algebra’s generator for a group SU{d) was used, where 
d > 2 is to define a set of operators responsible for transfer dynamics. For 
clarity, the construction procedure of SU{d) generators will be recalled - in 
the first step a set of projectors is defined: 


= \k){j\ = l<v,fj.<d. (31) 

The first suite of d{d — 1) operators from the group SU{d) is specified as 

0k j ^ pk,j pj,k^ ^k,j ^ (-32) 


and 1 < k < j < d. 

The remaining (d — 1) generators are defined in the following way 




r,r 



^pr+l,r+l 


(33) 


and 1 < r < (d — 1). Finally, the d^ — 1 operators belonging to the SU{d) 
group can be obtained. 
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